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ABSTRACT 

The Kepler mission results indicate that systems of tightly-packed inner plan¬ 
ets (STIPs) are present around of order 5% of FGK held stars (whose median 
age is ~ 5 Gyr). We propose that STIPs initially surrounded nearly all such 
stars and those observed are the final survivors of a process in which long-term 
metastability eventually ceases and the systems proceed to collisional consoli¬ 
dation or destruction, losing roughly equal fractions of systems every decade in 
time. In this context, we also propose that our Solar System initially contained 
additional large planets interior to the current orbit of Venus, which survived 
in a metastable dynamical configuration for 1-10% of the Solar System’s age. 
Long-term gravitational perturbations caused the system to orbit cross, leading 
to a cataclysmic event which left Mercury as the sole surviving relic. 

Subject headings: celestial mechanics — planetary systems — planets and satel¬ 
lites: dynamical evolution and stability — planets and satellites: formation 


1. Absent Planets 


Why aren’t there planets interior to Mercury? This question is particularly evident in 
light of the discovery of many multi-planet systems at distances <0.5 AU containing several 
Earth masses of material. Our answer is: there were, and Mercury is all that remains, which 
fits our Solar System into a framework where dynamical instability mercilessly consolidates 
or degrades close-in planets. 


The Kepler mission discovered many syst e ms of tighty-packed in ner planets (STIPs) fjFabrvckv et al. 


2014 iRowe et al.l 12014 iLissauer et al.l I20l4) . iLissauer et al.l (1201111 estimate that 
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Kepler stars host STIPs and Fressin et alj (120131 ) conclude that half of the Kepler stars have 
at least one 0.8 — 2 R 9 planet with orbital periods shorter than Mercury’s. The absence of 
such close-in planets in our Solar System is perhaps surprising; the surface mass density a 
profile for the minimum mass solar nebula (with radial dependence a oc a~ l or a -1,5 ) yields 
several Earth-masses of material inside 0.7 AU if the disk extends down to the ~ 0.05 AU 
distance where STI Ps are found and where the inner edge of gaseous protoplanetary disks 
are thought to be (Meyer et aL 1997 ). 


In contrast, Solar System terrestrial planet formation models require an inner edge 
to the pla netes imal disk at ~ 0.5 — 0.7 AU in order to reproduce Mercury’s small mass 


(jChambersll200ll; lHansenll2009l) and because the angular momentum of the t errestrial plane 


is inconsistent with accreting significant mass from interior to Venus’ orbit (Wetherill 19781). 
Histo rically this was not viewed as troubling, because chemical condensation modeling (jLewis 


19741 ) indicated that temperatures closer to the Sun would rise above the condensation 
temperature of any solids; these models’ high temperatures at 0.5 AU also seemed to explain 
Mercury’s metal-rich rich nature. 


However, ’dead zones’ close to stars may inhibit MHD turbulence (iLvra et al.l 2009j), 
reducing energy dissipation and tem peratu res in these opti cally thick regions; planet for¬ 
mation may sequester dust rapidly (Dzyurkevich et al. 201Q) resulting in the STIPs regions 
being undetectable in the protostellar SED. Other, ubiquitous disk pro cesses a lso promote 
the rapid formation of planetary building blocks very close to the star (Boley et ah 2014). 
Given a supply of solids interi or to 1 AU, accretion simulations show planets forming easily 
in these regions (jHansen & Murra y 201 31). Starting with the hypothesis that nearly all FGK 
stars form with a STIP, it is probable that such systems are dynamically metastable on a 
variety of timescales, allowing for planetary consolidation or destruction. In this context, 
our explanation for the Solar System’s current lack of large planets interior to ~ 0.7 AU 
is that our Solar System originally formed with a STIP at < 0.5 AU composed of a few, 
now absent, Earth-scale planets. Through multiple generations of catastrophic collisions and 
re-accumulations initiated by a dynamical instability between the original planets, we now 
have Mercury as the last remaining relic. 


2. Metastability in Planetary Systems 


Kepler host stars are ~1-10 Gyr old (iMarcv et al.l 12 0141 ); obviously STIP formation 
must allow long-term dynamical stability, even if s ome of the syst ems’ planets are nested at 


intervals barely beyond the stability requirements (ILissauer et al.l 1201 11). Some well-studied 


r 

i). 


systems are today on the edge of dynamical instability (ILissauer et al.l 2013; Deck et ah 
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2012; Lissauer et ah 2012). This seems completely reasonable: planet formation gradually 
combines dynamically unstable protoplanets, so evolving systems will rarely transition from 
being ‘highly coupled’ to ‘extremely oversta ble’. Thus planetary systems should naturally al¬ 
ways be in a state of metastability f Laskar 1996h . Our Solar System itself is only metastable; 
the terrestrial planets’ orbits are chaotic ( Laskarl 1989), and Mercury has a 1 % chanc e of cre¬ 
ating large-scale chaos on 5 Gyr timescales (lLaskarlll996l: lLaskar fe Ga.stineaul 120091 ) . STIPs 
should exhibit similar metastability, with many systems metastable on the lifetime of their 
star, while others reached orbit crossing in the past. In this framework, the STIPs frequency 
found by Kepler represents a lower limit on their formation probability because we only see 
the Gyr stable systems. The absence of STIPs around many stars could be due to the earlier 
collapse of a metastable planetary arrangement, leaving one or no detectable planets at short 
periods. 

To explore metastable states in STIPs, we preformed a large suite of numerical integra¬ 
tions based on the observed, presumably Gyr stable, Kepler STIPs. We generated systems 
with architectures similar to the known systems by using the observed semimajor axes and 
plane tary radii, calculating planetary masses using the relationship M p ~ M E (R p /R E ) 2m 


(iLissauer et al. 2011). We randomized the initial orbital angles, assumed nearly coplanar 
orbits (mutual i < 1.5°), and assigned random initial eccentricities e 0 =0-0.05; if the ob¬ 
served systems had measured maximal e, we assigned eo =0-e max ■ The real Kepler systems’ 
eccentricities are weakly observati onall y constrained (if at all), but the range we consider is 
consistent with the observations (Fabrycky et al. 2014). We found that the range of eo for 
the architectures we explored was unimportant; instability probability is not correlated with 
eo (within our chosen range). 

We integrated analogs of 13 observed Kepler systems with more than 4 planets: Kepler- 
102, Kepler-84, Kepler-90 (Lissauer et al. 2014); Kepler-107, Kepl er-169, Kepl er-292, Kepler- 

20141): Kepler-11 (iLissauer et al. 20131): Kepler-62 ( Borucki et al. 


223, Kepler-26 flRowe et al. 


2013 ): Kepler-85 ( Ming et al.ll2013 ): Kepler-20 ( Gautier et al.ll2012 ): and Ke pler-33 (Lissauer et al. 


2012). These analogs were evolved using the Mercur y orb ital int egrator (IChambers 19991 ) 
with an approximate inclusion of general relativity (ILissauer et al.l 2 011 ). We integrated 
20 analogs of each Kepler STIP for 10 Myr or until a physical collision occurred between 
two planets. For Kepler-169, 292, and 84, no analogs had a collision within 10 Myr from 
our initial conditions; additional 10-Myr simulations with eo =0.05-0.1 also produced no 
collisions, indicating that these system architectures are likely stable on long timescales. In 
our framework these may represent systems which have had a planetary consolidation in 
the past, or which are simply unstable on longer, 1-10 Gyr timescales. This does not imply 
that the other ten observed systems are currently unstable on 10 Myr timescales. We are not 
investigating the stability of the exact observed systems, which have been shown to be stable 
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for ~10 8 yr assuming initially circular orbits flLi ssauer et al.l 2 01 ll; Fabrycky et ah 2014). We 
are instead interested in the range of possible behavior for STIPs with planet masses and 
spacings similar to the observed STIPs. 

For the ten systems that had collisions, we ran large suites of similar integrations for 
100-150 Myr to explore the questions: (1) how do systems evolve as instability is approached? 
and (2) how are the instability timescales distributed? Fig. |T] shows typical evolutions for 
analogs which exhibited a first planetary collision after 45-140 Myr of metastability. Out of 
the ~600 system analogs we integrated, half resulted in a collision between planets. What 
may be surprising is that there is no obvious sign of coming instability; as Fig. [Q shows, 
the systems maintain e ~ eo for nearly the entire duration before suddenly transitioning 
to orbit crossing. Initial planetary eccentricities have no systematic effect on when the 
instability occurs, but even tiny differences in initial conditions can produce vastly different 
instability timescales. Pot ential sources of chaos in exoplanet systems have recently been 
discussed by iQuillenl (1201111 and iDeck et al.l (1201311 ; how sys tems can evol v e to a suddenly 
more chaoti c state at an unpredi ctable time are discussed bv iLaskarl (1199611 . iLithwick fc Wu 
( 2014 1. and Batygin et al. ( 2015 1. 


Our experiments show that instability timescales in these systems are distributed such 
that equal fractions of the systems go unstable (reach a first planetary collision) in each 
decade in time (Fig. [2]). This logarithmic decay is not unknown in dynamical systems 


eg., 


Holman & Wisdom (19930) and is presumably related to chaotic diffusion. After a 
brief, relatively stable initial period, the systems hit instability at a rate of ~20% per time 
decade, with half of the systems still intact at ~100 Myr. The exact decay rate may be 
influenced by our usage of the current Kepler STIPs (perhaps the most stable); however if 
this decay rate held, at ~5 Gyr 5-10% of STIPs would not yet have reached an instability, 
in rough agreement with the observed STIPs frequency. 

Mercury being a remnant of a p reviously unstable system fits nicely with its current 
instability timescale (ILithwick &; Wul 12014); it would also n ot be the first suggestion of 
a metastable Solar System configuration (Gomes et al. 120050 . If our Solar System once 
contained a metastable set of planets interior to Venus, one can eliminate the artificial 
inner disk ed ge used in terrestrial p lane t formation models; the difficulty in making Mercury 


analogs (IChambersI l200ll; lHansen 


20091) would then be explained by the fact that Mercury 


is a collisional remnant. Test integrations determined that the orbital evolution of the three 
outer terrestrial planets (Venus, Earth and Mars) are unaffected on 500 Myr timescales 
by the presence of four additional planets totaling 4M 0 in mass with a < 0.5 AU. We also 
observe in our STIP simulations that when inner planets experience instability, the outermost 
planet (Venus analogs at > 0.5 AU) is often unperturbed. Thus it is not unreasonable that 
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Kepler-11 analog 



Fig. 1.— Evolution of semimajor axis, periapse distance, and apoapse distance for fonr STIPs 
analogs. Sometimes the innermost two planets collide leaving the outer planets relatively 
undisturbed (top panels) while other systems show close encounters between three planets 
(bottom panels). 
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fractional time 


Fig. 2.— Fraction of STIP analogs that have not experienced a collision as a function of 
integration time. The dashed line shows a slope of 20% loss per decade in time. 

Venus, Earth, and Mars could escape large-scale orbital effects as the Solar System’s STIP 
disintegrates. 


3. Consolidation and/or Destruction 

Once instability is initiated, the possible end-states of STIPs will fall along a continuum 
with two extremes: 1) consolidation , where almost all of the initial mass ends up in a smaller 
number of planets, or 2) destruction, when <10% of the STIP’s mass survives. We propose 
that our Solar System reached the destructive end state, but other systems consolidated an 
initial many planet STIP into fewer, more massive short-period planets. The destructive 
end state is made possible by the extreme collision speeds which can occur for such close- 
in orbits; further from the star, the ratio of typical impact speeds, v imp , to mutual escape 
speeds, v esc = y/2G(M 1 + M 2 )/(R\ + i?/), are low enough that accretional/consolidational 
processes are more likely. We note that most of the literature on planet-scale collisions has 
understandably focused on the context of accretion near 1 AU; at ~0.1 AU impact speeds 
rise by factors of ~3, greatly increasing the likelihood of erosive collisions. 

Fig. [3] shows an example of how a system might evolve after an instability. Like in 
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Fig. IU there is initially no macroscopic evidence of instability. In the last 10% of that phase, 
a transition occurs where several planets begin to interact, leading planets b and c to collide. 
We had (incorrectly) postulated that essentially all analogs would reach their first collision 
as a result of an instability between only one pair of planets; we find instead that ~ 25% 
of the analogs show close encounters between 3 or more planets before the first pairwise 
collision occurs. 


We continued some simulations beyond the first collision to study generic features of 
the subsequent evolution. To do this, we assumed perfect, inelastic merging of the colliding 
planets (which is unlikely to be a good approximation for reasons discussed below). The 
outcomes subsequent to the first collision are highly ergodic; in Fig [3J, the 6 planet system 
consolidates to 3 planets which then remain stable for at least 100 Myr. 


Marcus et al.1 (1200911 showed that the ratio (v imp /v esc ) is important in determining the 


frequency of erosive collisions, which are more common when (vi mp /v esc ) > 2. Similarly, 


Stewart fe Leinhardti (2012) showed collision outcomes as a function of ( Vi mp /v esc ) and im¬ 


pact ang le broken down into accretion, catastrophic collision, and ‘hit and run collisions’ 
( Asphaug et all 20061) in which grazing collisions liberate some mass. Even outside the 
catastrophic disruption regime, these authors point out that significant fractions (~10-15%) 
of impacting mass in any collision is likely dispersed into small debris. 


We compiled ( v irnp /v esc ) values from our simulations (Fig. [I]); for systems initially spaced 
by 10-30 mutual Hill radii (like the Kepler analogs), if two comparable-mass neighboring 
planets at 0.1 AU are excited to mutual crossing, then e ~ 0.2 at the time of first crossing, 
yielding encounter speeds v enc ~ e Vkepier ~ 20 km/s. Super-Earths with ~ 1.6 Rq have es¬ 
cape speeds near this value, so when the first pair of planets collide Vi mp = a / v\ nc + v\ sc ~ 1, 
as Fig.[]]confirms. These first collisions are likely consolidational. If this was the whole story, 
one might think that STIPs gradually pairwise combine in a dominantly accretional environ¬ 
ment. However, our simulations show that subsequent collisions often occur at much higher 
impact-speed ratios (Fig. [4]). The previous excitation in the system results in higher speed 
second collision, either nearly immediately (vertical lines in Fig. [4j) or after a metastable 
delay (diagonal lines). Many of these second collisions rise into the more erosive regime; we 
propose multi-collision excitation is how some systems enter the destructive regime. 


4. Mass loss mechanisms 

After instability, many STIPs may experience only a sequence of low-speed, accretional 
impacts between nearest-neighbor planets that produce small amounts of rapidly-eliminated 
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Fig. 3.— Instability in a Kepler-11 analog. An initial period of apparent stability (upper 
panel) transitions to a perturbed state with three collisional events (modeled as perfect con¬ 
solidations) before reaching a final, 3-planet metastable state (middle panel); the shaded 
period between the first and second collisions involves four planets in dramatic ‘orbit swap¬ 
ping’ (bottom panel). 
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Fig. 4.— Ratio of collision speeds to mutual escape speed vs time for ST1P analogs. First 
collisions (black dots) are dominantly accretional {v„ np /v esc < 2). The second collision in a 
given system (red dots connected to the first event by a line) clearly shows that sequences 
of collisions allow more erosive conditions. The highest collision speed was 75 krn/s 
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debris. Such systems consolidate down to fewer planets, concentrating the mass spectrum 
toward larger planets. Indeed the distribution of planet sizes in the Kepler systems with 
3 or more planets show a trend toward larger planets at lower multiplicity (not formally 
significant). Even in consolidational systems, instabilities can produce moderate inclinations, 


causing some surviving planets to become undetectable in transit surveys (IJohansen et al. 


2012 ). 


Our conclusion that some fraction of STIPs reach a first instability and then consolidate 
or degrade is independent of the hypotheses that follow regarding possible significant mass 
loss. However, some of our STIPs analogs experience heavy perturbation; such systems are 
candidates for substantial mass loss either through continuous bursts of debris dispersal or 
even large-scale planetary elimination. We divide this into four size scales: 


1. Dust below the blow-out limit (~1-10 microns) hyperbolically leaves the STIP region 
in just months. This mass-loss mechanism is very efficient if planetary-scale events 
directly generate large amounts of dust or if dust is produced by a cascade of smaller 
collisions in the aftermath of each major event. 


PR drag can cause cm-scale particles to spiral from 0.1 AU down to the star on 
timescales of 100 kyr. However, if a collisional cascade produces a considerable amount 
of small d ebris, then t he tim escale for debris self-collision is shorter than the PR drag 
timescale fjMelis et al, 2012) and particles cannot inspiral before being reduced to dust 
and blown out. 


The above two processes might be inefficient if ~0.1 M e of cm or smaller debris is produced 
in any singl e ev ent because the optical depth to the star exceeds unity and the disk could 
self-shield (Gladman & Coffey 2009), shutting down radiation effects. A competition can 


occur between the timescale for the largest remnant to sweep up debris and the timescale for 
debris to self-collicle and grind down to the PR and dust blow-out scales. Moderate events 
(~ 100 km scale) in that collisi onal cascade produ c e sudden spikes in dust that quickly decay, 
perhaps like those observed bv lMeng et al.l (120141) : ISong et al.l (120051) . Large-scale planetary 
violence finishes within ~ 5 x 10 4 yr (see Fig. [3]), so only 10 -5 of mature field stars would 
show these sudden bursts of hot-dust excess. Because of the ‘equal fraction per decade’ 
instability behavior, samples of younger stars would have hot-dust probabilities inversely 
proportional to their age. 


3. For meter to 100 km objects, evolution is largely driven by repeated gravitational 
scatterings by the largest remnant (s). With high relative orbital speeds, gravitational 
focussing (which enhances re-accretion) is minimized. High results in erosive 
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impacts for most impact angles ([Marcus et al. 2009), potentially removing, rather than 
adding, mass from the remnants. The 50-100 km/s impact speeds occurring this close 
to the star may enhance vapor production, hindering ejecta retention relative to slower 
impacts out near 1 AU. In the Solar System, debris interior to ~ 0.5 AU has very short 
collisional lifetimes and is subject to removal via Yarkovsky drift , con sistent with the 
current lack of km-scale and larger debris in this region (tVokronhlickv et ah 20 00 ). 


4. Close to a star, pure ejection is unlikely because v esc /vk ep ier is so small. However, 
secular interactions in a post-instability system could push planets or debris to star¬ 
grazing (e ~ 1) on Myr timescales. Although stellar impacts are unobserved in our 
simulati ons, e pumping via secula r resonances is a known p henomenon in our Solar 
System I Na/monni fe Murray 1999 ; Laskar fe Gastineau 2009h and often relies on the 
presence of exterior giant planets which are unseen (and thus unmodeled) in the Ke¬ 
pler systems. While the Kepler planets are obviously not today near secular resonances, 
post-instability evolution could change this; Fig. |3h illustrates scattering planets ex¬ 
ploring a large range of a, e, i space, a near-perfect algorithm to find secular resonances. 
It is difficult to assess the probability of significant secular evolution, but we ran sim¬ 
ulations of pairs of ~Earth-mass planets evolving in our current Solar System (minus 
Mercury) on low-e orbits form 0.1 — 0.4 AU; this configuration is a plausible outcome 
of a recently consolidated 3-planet STIP. We find many configurations where these 
planets’ eccentricities grow to e = 0.4 — 0.9 on ~ 1 — 10 Myr timescales, with the inner 
planet’s perihelion distance sometimes dropping to just a few Solar radii. Although 
achieving e ~ 1 is rare, moderate secular eccentricity growth can promote subsequent 
high-speed collisions between the remaining planets. 


We thus postulate that our Solar System originally had a STIP of 3 or more now-absent 
planets totaling a few Earth masses. An instability initiat ed a s equence of collisions (as 
opposed to a single collision of a ~ 0.2 M e body, Ben z et al. ( 19881B . which allowed the sys¬ 
tem’s excitation to the destructive regime; suc h a process concentrates iron into the surviving 
remn ants, explaining Mercury’s high density (jStewart fc Leinhardtll2012t lAsphaug fc Renter 
2014 b Mercury’s current e ~ i (radians) ~ 0.2 Q and Mercury’s current ‘survivor’ metasta¬ 
bility timescale of ~5 Gyr would naturally result from this scenario. Additionally, the 


transplant of iron-meteorite parent bodies to the asteroid belt (IBottke et ah 2006) from the 
< 0.5 AU region in this scenario can easily acco mm odate their rapid initial accretion and 
evidence for grazing protoplanetary impacts (Goldstein et ah 2009). 


1 That is, e 


Vesc/Vkep ~ (20 km/s)/(100 km/s) 
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Much work beyond this Letter is required to explore the details of this general frame¬ 
work. For example, the probability of destructive end states should be consistent with the 
de-biased estimate that half of mature FGK stars have no visible planet with R > 1 R® and 
P < 88 days ( Fressin et ah 2013 b 


5. Mercury and the Lunar Cataclysm 


The decadal nature of the instability makes it plausible that the penultimate metastable 
state lasted ~ 0.5 Gyr (one tenth Mercury’s current metastablility timescale), allowing the 
additional speculation that this last instability, 4 Gyr ago, was responsible for the “Lunar 
Cataclysm” (reviewed by H ar t mann et al. (120001 )). An instability transitioning into rapid 
pl anetary des truction (in ~1 Myr) would spread debris throughout the inner Solar System. 


Glad man fe Coffevl (120091 ) estimated that 10-20% of m to 100 km debris originating near 
current Mercury would strike Venus, with 1-4% impacting Earth (~0.1% strikes the Moon). 
The Earth’s impact rate would peak ~1—10 Myr after the event and decay on ~30 Myr 
timescales as Mercury and Ven us ab sorb most of the debris; this is a plausible match for the 
cataclysm’s final stages (ICuk et al.l I2010T) . A bottom-heavy size distribution for the 1-100 
km debris could explain the recent finding (IMinton et ah 20151 ) that a main-belt asteroid 
source would produce too many impact basins during the cataclysm. ST1P debris would 
likely be mostly silicate-rich mantle material similar but not identical to main-belt asteroid 
compositi ons, consistent with cataclysm impactor compositions inferred via cosmochemical 
means (JJoy e t ah 2012). The smallest dust (being blown out hyperbolically) could impact 
the Earth-Moon system. We estimate that 10 -11 of the departing dust would strike the 
Moon, at v imp ~30 km/s. If any dust or meteoroid projectiles were retained, fragments 
might be found in regolith breccias compacted during the cataclysm epoch. Compared to 
traditional cataclysm hypotheses, this scenario yields significantly higher impact rates onto 
Venus, with potentially significant implications for its evolution. 
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